Neutrino Pair Emission from Hot Nuclei During Stellar Collapse 
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We present shell-model calculations showing that residual interaction-induced configuration mix- 
ing enhances the rate of neutral current de-excitation of thermally excited nuclei into neutrino- 
antineutrino pairs. Though our calculations reinforce the conclusions of previous studies that this 
process is the dominant source of neutrino pairs near the onset of neutrino trapping during stellar 
collapse, our shell-model result has the effect of increasing the energy of these pairs, possibly altering 
their role in entropy transport in supernovae. 

PACS numbers: 



I. INTRODUCTION 

In this paper we study the role of the nucleon-nucleon 
interaction in the process of de-excitation of hot, excited 
nuclei into virtual Z 0, s and neutrino-antineutrino pairs. 
This process is likely the dominant source of neutrino 
pairs in collapsing stellar cores [ij-y]. The energy of the 
neutrinos in these pairs, set in part by nuclear structure 
considerations, can be an important determinant of en- 
tropy transport in core collapse supernovae. In turn, the 
entropy figures prominently in the nuclear composition, 
neutronization history, and initial shock energy in super- 
nova models. 

The entropy per baryon in the collapsing stellar core is 
low and, as a result, most nucleons reside in large nuclei 
and there are very few free protons [3, @ ■ The paucity of 
free protons has the effect of suppressing the overall elec- 
tron capture rate @ , yielding a greater electron fraction 
Y e (electrons per baryon). Most pressure support within 
the core comes from electron degeneracy, so higher Y e 
during collapse implies a larger pressure and, hence, a 
larger homologous core mass. The mass of this inner core 
determines the initial energy of the post-bounce shock: 
a more massive inner core yields a stronger initial shock. 
The strength of the shock, the mass of the core above 
the shock, and photo-dissociation of heavy nuclei in this 
outer core (all determined in part by Y e ) are important 
parameters in the supernova explosion process [6l4l8l| . 

In the epoch near neutrino trapping, when the 
core density is ~ 10 12 gcm -3 , the electron fraction is 
Y e 0.32 [(J, giving an electron Fermi energy \i e « 

51.5 MeV (pi2 Y e ) 1/a w 35MeV, where the density is 
scaled as p%2 = p/10 12 gem -3 . The temperature of the 
core is in the neighborhood of T ~ 1 MeV to 2 MeV, so 
the electrons are highly degenerate. (In this paper we 
use natural units and set h = fcb = c = 1.) Energy 
emission from the core via neutrinos helps to maintain 
low entropy, but at a core density of p\2 ~ 1, high en- 
ergy neutrinos are trapped by neutral current coherent 
scattering on nuclei. However, since the cross-section for 
this process varies as the square of neutrino energy, low 
energy (E u < 10 MeV) neutrinos may escape, carrying 



away entropy and possibly, depending on the process, 
lepton number. 

There are a number of ways to produce low energy neu- 
trinos in the core: inelastic down-scatter of neutrinos on 
electrons [lj| [2(|; electron neutrino-pair bremsstrahlung; 
plasmon decay, etc. (21- 29] . However, since these pro- 
cesses involve electron interactions, most of them are sup- 
pressed by the extreme electron degeneracy. In these con- 
ditions there simply isn't much phase space available for 
electron down-scatter. 

This leads us to explore low-energy neutrino produc- 
tion mechanisms that do not involve electrons. Simple 
neutral current neutrino-nucleon down-scattering tends 
to be ineffective in this regard. This is because this pro- 
cess is roughly conservative: the nucleon mass is large 
compared to the typical electron capture-generated neu- 
trino energies, E u ~ 25 MeV. Free nucleon neutrino-pair 
bremsstrahlung (3fj| - l32l |. a key process in neutron star 
cooling, is less effective here because of phase space con- 
siderations and because there are few free nucleons in 
the low entropy conditions that favor large nuclei during 
stellar collapse. 

Indeed, there are analogs of these two processes for nu- 
cleons that reside in the large nuclei characterizing the 
neutrino trapping epoch, and these are not subject to the 
limitations of their free nucleon cousins. These analogs 
are: inelastic down-scattering of energetic neutrinos on 
nuclei; and neutrino pair emission from thermally-excited 
nuclei 0, [H, [33| . The first of these processes has a rela- 
tively large cross section, as this channel is what Ref. [l| 
termed an "up" transition, where the nucleus acquires 
energy from the neutrino, implying that nucleons tran- 
sition to higher nuclear excited states, where they are 
relatively less Pauli blocked. These processes can have 
important implications for supernovae [l], l34j - [37| ] . 

In contrast, the second of these, the nuclear de- 
excitation into neutrino pairs channel shown in Fig. [TJ 
is a "down" transition, subject to more nuclear Pauli 
blocking, and therefore possessing considerably less nu- 
clear weak interaction strength on average than the neu- 
trino inelastic down-scatter channel. Nevertheless, the 
de-excitation process has some unique features: in princi- 
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pie it may produce lower energy neutrinos than the down- 
scattering channel and, should these escape the core, en- 
tropy will be lost but electron lepton number will not be. 
As shown in Ref. Q, hot nuclei can also de-excite into 
neutrino pairs through a virtual plasmon (photon col- 
lective mode in the plasma), and this process has been 
argued to lead to large enhancement factors in nuclear 
neutrino pair emission in the first forbidden channel. As 
we will show, our nuclear structure considerations also 
impact this channel. All of these issues depend to some 
extent on the nuclear physics of down transitions, and so 
this is where we concentrate in this work. 




A* 



FIG. 1: Neutral current neutrino pair emission from an ex- 
cited nucleus A*. 




FIG. 2: Thermally populated nuclear state with excitation 
energy E, de-excites via virtual Z emission to a final state 
with excitation energy E/. 



In Section II we discuss the nuclear and phase space as- 
pects of de-excitation into neutrino pairs. Nuclear shell- 
model considerations are discussed in section III. In Sec- 
tion IV we discuss results, and in Section V we give con- 
clusions. 



II. DE-EXCITATION RATES AND NUCLEAR 
STRUCTURE 

A. Large, Highly-Excited Nuclei During Core 
Collapse 

Strong and electromagnetic interactions are fast 
enough that the material in stellar collapse can be in 
thermal and chemical equilibrium, i.e., Nuclear Statisti- 
cal Equilibrium (NSE) . The weak interaction also is driv- 
ing toward equilibrium (beta equilibrium) at the epoch 
of neutrino trapping, but has not yet arrived there. The 
pioneering work by Bethe et al, Ref. showed that the 
entropy-per-baryon in a collapsing iron core is s s=s 1 and, 
as a consequence, most all nucleons will reside in large 
nuclei. Minimizing the free energy for typical conditions, 
for example with p 12 ~ 1 and T = 1 MeV to 2MeV, 
yields a mean nuclear mass A ~ 100. 

The mix of nuclei in NSE in these conditions is exotic. 
These huge nuclei will have fair neutron excess, because 
Y e < 0.4, implying neutron-to-proton ratios n/p > 1.5. 
Moreover, because the nuclear level density is high, these 
nuclei will be sitting at high excitation energies. To see 
this we can treat the nucleons in the nucleus as a Fermi 
gas. Using the familiar Bethe approximation [38[ for the 
nuclear level density,an estimate of the average nuclear 
excitation energy is [J] 



(E) w aT 2 



(1) 



where the level density parameter is a « A/ 8 MeV 1 . 
For example, with A = 100 and T = 2 MeV, Eq. [T] im- 
plies an excitation energy (E) ~ 50 MeV. The expression 
in Eq.[T]is easily understood: The number of nucleons ex- 
cited above the nuclear neutron and proton Fermi levels 
will be r; aT, and in thermal equilibrium each nucleon 
so excited will have an average excitation ~ T. 



B. Nuclear De-Excitation Rates 

Consider an excited nucleus dropping down to a lower 
excitation energy via virtual Z° emission, as depicted in 
Fig.O This is similar to a nuclear Ml gamma transition. 
The de-excitation rate [fi[3^ from initial nuclear state \i), 
with excitation energy Ei, to final nuclear state |/), with 
excitation energy Ef, is 



A 



if 



^lA {A E) 5 B(GT) lf 



60 7T 

1.71 x 10~ 4 



/ AE \ 



MeV J B ^ 



(2) 
(3) 



where Gf is the Fermi constant, <?a ~ 1-26 is the ax- 
ial vector coupling constant, AE = \Ei — Ef\ (hereafter 
referred to as the transition energy) is the difference be- 
tween the initial and final state nuclear excitation en- 
ergies, and B{GT) l} = |(/||S fc (^ t z ) k \\i)\ 2 / + 1) is 
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the reduced transition probability associted with the ax- 
ial vector operator. The matrix element connects initial 
nuclear state \i) with final nuclear state |/}. Here ~a is 
the Pauli operator and t z is the z-componcnt of nuclear 
isospin. For the nuclei we consider here, only the axial 
vector matrix element is significant: when we neglect the 
relatively small Coulomb effects, the nuclear part of the 
Hamiltonian commutes with isospin operators {e.g., T z ), 
and the weak vector matrix element |(/|T z |i)| is zero. 

The corresponding neutrino-plus-antineutrino energy 
emission rate, A*/, for this transition is the product of the 
de-excitation rate and the transition energy. Whether 
or not the neutrinos carrying this energy escape from 
the star without scattering, thereby turning the energy 
emission rate into an energy loss rate, depends on many 
factors, most especially the neutrino energies. 

The total energy emission rate for an excited nucleus in 
initial state \i) is the sum of the energy emission rates to 
all accessible final states, and can be viewed as a function 
o{E h 

K ot (Ei)= J2 \Ef-Ei\X if . (4) 

/, E f <Ei 

The total overall energy emission rate for the entire nu- 
cleus follows on performing a population index-weighted 
sum over all initial states i, 

1 r°° 

A tot = J2 ^ A* ot ~ ^ / P (Ei) e- E >/ T A* ot dE h (5) 

where P, ; = (2 J, + 1) exp(-^/T)/Z is the popu- 
lation index for state i, with Jj the spin of level i, 
Z = J^i (2 Ji + 1) exp(— Ei/T) is the nuclear partition 
function, and p(Ei) is the nuclear level density at excita- 
tion energy Ei. 

For the thermodynamic conditions relevant for NSE 
near p\2 ~ 1, we can get a crude estimate of A to t by 
simply evaluating A' ot at the mean excitation energy for 
temperature T, i.e., taking Ei — (E), 



A tot «Ar = (6) 

The rationale for this approximation is that while the 
nuclear level density rises nearly exponentially with exci- 
tation energy, the Boltzmann factor in Eq. [5] falls ex- 
ponentially with this energy, so that their product is 
strongly peaked at (E). 

However, using this rough approximation is problem- 
atic. The level density is high near (E), and there will 
be many different kinds of nuclear many-body states with 
e.g., different spins and isospins, but all with roughly this 
excitation energy. Therefore, choosing a single represen- 
tative state is not possible. 



III. SHELL-MODEL CONSIDERATIONS 
A. Approaches to the Problem 

Evaluating the energy emission rates in Eq.|4]and Eq.[5] 
for a nucleus with nuclear mass number A ~ 100 at a 
mean excitation energy ~ 50 MeV is clearly impractical 
with conventional nuclear structure techniques tailored 
to capture low excitation energy physics. The shear size 
of the problem, some two dozen particles excited above 
the Fermi surface in a mass ~ 100 nucleus with all of the 
fp-, gd-, and gh-shells in play, precludes this route. 

There are two possible alternative approaches: (1) 
Treat the ~ aT nucleons excited above the Fermi sea 
as nearly free particles within a dense environment, with 
appropriate phase space modifications, and then calcu- 
late the neutrino-pair bremsstrahlung rates for these; and 
(2) Exploit the fact that each nucleon excited above the 
nuclear Fermi level has only a relatively small amount 
of energy (~ T) , so that conventional shell-model treat- 
ments are efficacious, at least for nuclei with low enough 
mass that the problem is tractable computationally. The 
first of these, by treating valence nucleons as plane 
waves @, will tend to overestimate [l[ the nuclear weak 
strength available, but has the advantage that it would go 
smoothly to the homogeneous matter limit when nuclei 
merge at high density (pi2 > 10). 

Here we will take up the second approach, in part be- 
cause it has the advantage of getting a better handle on 
the weak nuclear strength, nuclear structure effects, and 
energetics. The latter is an especially critical issue since 
the neutrino energy emission rates scale like six powers 
of the transition energies. Therefore, ascertaining how 
e.g., configuration mixing and particle-hole repulsion act 
is important. However, we will have to model nuclei 
with lower masses, generally sd-shell and fp-shell species, 
rather than the mass ^100 nuclei that neutrino trapping 
NSE conditions pick out. At best, this approach will al- 
low us to see trends that may at some point facilitate 
extrapolation of these considerations to enable estimates 
for rates from heavier nuclei. 



B. Extension to Heavy Nuclei 

Even calculating the de-excitation rate from a single 
level in a "small" nucleus is a challenging and unusual 
nuclear structure problem, in part because matrix ele- 
ments between highly excited states are required. We 
have approached this using a conventional nuclear shell 
model with the usual filled closed core of nucleons in 
low-lying single-particle states plus valence nucleons in a 
model space. We then employ the usual Lanczos itera- 
tion with an appropriate nuclear Hamiltonian, with the 
twist that our start vector is taken to be a very highly 
excited state. 

Using the shell-model code Oxbash [39] , we performed 
a full sd-shell calculation of 28 Si using the USDB Hamil- 
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tonian 40] (closed 16 core with 12 valence nucleons in 
the Id and 2s shells). We performed a full fp-shell calcu- 
lation of 47 Ti using the GPFX1 Hamiltonian [Zl[ (closed 
40 Ca core with 7 valence nucleons in the If and 2p shells). 
Finally, we performed a truncated fp-shell calculation for 
56 Fe using the GPFX1 Hamiltonian, only allowing up to 
2 valence protons and up to a total of 4 valence nucleons 
to occupy single-particle states above the zero-order (no 
residual interaction) ground state configuration. Some of 
the sd calculations were carried out with the NuShellX 
code H|. 

When experimental Gamow- Teller beta-decay 
strengths are compared to the results obtained from 
calculations in the sd and pf model spaces it is observed 
that experimental strengths are uniformly reduced 
relative to theory by a factor of 0.5-0.6 [U, 51]. This 
"quenching" is mainly due to second-order configuration 
mixing induced by the short-ranged part of the tensor 
interaction {451 [iq ] . We assume that the same quenching 
applies to transitions from the excited states and we use 
a reduction factor of r = 0.50. 

We examined transition strengths and energy loss rates 
over a range of excitation energies from to 40 MeV 
in 28 Si and 56 Fe and at 23 and 27.6 MeV excitation in 
47 Ti. The strength distributions for all three nuclei at 
27.6 MeV excitation are shown in Fig. [3] along with the 
distribution obtained by averaging the strength as a func- 
tion of transition energy over all three nuclei. While the 
details of the shapes of the strength distributions vary be- 
tween nuclei, the essential feature of a central peak with 
a long tail out to transition energies of 15 or 20 MeV is 
consistent. 



0.18 
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Transition energy E, - E, (MeV) 

FIG. 3: Neutrino-pair emission transition strength for 28 Si, 
47 Ti, and 56 Fe from initial states at 27.6 MeV excitation 
shown as functions of the difference between initial excitation 
energy Ej and final excitation energy E/. 

To obtain energy emission rates per nucleon as func- 
tions of excitation energy, shown in Fig. 21 we applied 
Eq. 0] to each nucleus, then divided by A. With an eye 
toward extension to large nuclei, the key observation is 



that the energy loss rate per nucleon is strongly depen- 
dent on excitation energy, but nearly independent of nu- 
cleus, despite the considerable differences in the models 
used for each nucleus. 




5 10 15 20 25 30 35 40 
Excitation energy (MeV) 



FIG. 4: Energy loss by neutrino-pair emission for Si, Ti, 
and Fe shown as functions of excitation energy. 

To find temperature as a function of mean excitation 
energy and nucleus, we invert Eq.[TJ This, along with the 
approximation in Eq. [6] gives emission rate per nucleon 
as a function of temperature, shown in Fig.[5j Also shown 
is the result for 28 Si, but with the temperature computed 
from excitation energy as though it had the same mass 
number as 56 Fe. As can be seen from Eq.[T] this amounts 

1/2 

to scaling the temperature of 28 Si by a factor of (||) 
This allows us to compare the 28 Si (our most realistic 
model) and 56 Fe (our most astrophysically relevant nu- 
cleus) results directly as functions of temperature: the 
comparison of the scaled results in Fig. [5] is equivalent to 
what is shown in Fig. |4l We use this method of tempera- 
ture scaling extensively throughout this paper; it will be 
indicated in each case. 



IV. RESULTS 

We found that the weak strength in transitions from 
highly excited initial states is spread significantly in en- 
ergy. The 28 Si results for the distribution of axial vector 
strength (B(GT): squared matrix element) with transi- 
tion energy and the corresponding neutrino-pair energy 
spectrum are shown in the lower and upper panels, re- 
spectively, of Fig.[SJ This plot shows both "up" transition 
strength, corresponding to positive values of Ef — Ei, and 
"down" strength with negative values of E / — E{ appro- 
priate for de-excitation into neutrino pairs. At 20 and 30 
MeV excitation, there is obviously more strength in the 
up than in the down channel. However, the up-strcngth 
distributions for both 30 and 40 MeV are not complete 
(artificially cut off) due to a truncation in the calculation 
at 10,000 Lanczos iterations. We obtained these distri- 
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FIG. 5: Energy loss by neutrino-pair emission for Si, Ti, 
and 56 Fe shown as functions of temperature. The dashed line 
shows the 28 Si result with the temperature scaled by a factor 

of (ff) 1 ^ 2 , allowing the 28 Si and 56 Fe results to be compared 
directly at equivalent excitation energies, as in Fig. [4] 



butions by averaging over ten states with Jj = 5 near 
each indicated initial excitation energy. Ji = 5 was cho- 
sen because the sd shell (2 J + 1) state density peaks at 
this spin value. 

According to Eqs. [3] and [U the strength distribution 
as a function of transition energy can be multiplied by 
six powers of the transition energy and a constant to give 
the contribution of a given transition energy to the overall 
neutrino-pair energy emission rate, i.e., the neutrino-pair 
energy spectrum. Although the actual strength distribu- 
tion is skewed toward lower energy transitions, weight- 
ing with six powers of transition energy clearly favors 
larger transitions. Summing over transition energy (e.g., 
Eq. @| gives the total neutrino-pair energy emission rate 
per baryon. Table Q] shows total energy loss rates per 
baryon computed in this fashion for several nuclei over a 
range of excitation energies and temperatures (rescaled 
as in Fig. [SJ. 

Our calculations show that the central peaks in the 
Fig. [6] strength distributions stem primarily from lateral 
(no spin flip) transitions. Such transitions do not change 
the single- particle energy of the transitioning nucleon. 
The wings of these strength spectra at larger transition 
energy come mostly from nucleon spin-flip transitions. 
Because configurations that result from a spin flip have a 
lower zero-order energy as a result of particle-hole repul- 
sion and spin-orbit splitting, they tend to be more readily 
mixed down to lower excitation energy than their coun- 
terparts stemming from no-spin-flip, lateral transitions. 
This became abundantly clear when we computed 28 Si 
with no spin-orbit splitting in the nuclear Hamiltonian, 
treating the i + 1/2 and i — 1/2 single-particle states as 
though they have the same energy. The results of this 
computation are shown in Fig. where it is readily seen 
that all of the strength is concentrated in a central peak, 




-30 -20 -10 10 20 

AE (MeV) 

FIG. 6: Bottom: Transition strengths for 28 Si from states at 
20, 30, and 40 MeV as functions of transition energy. Top: 
Energy emission rates per nucleus via neutrino pairs for 28 Si 
as functions of transition energy, i.e., spectra of emitted neu- 
trino pairs. 



drastically reducing the rate of neutrino pair production. 
Since the spin-orbit splitting is due to the nuclear sur- 
face, we conclude that the Gamow- Teller down-strength 
is greatly enhanced when the baryons are confined to nu- 
clei. 



V. DISCUSSION AND CONCLUSIONS 

We have found that transitions between spin-orbit 
partners account for the bulk of the spread to lower exci- 
tation energies (hence, larger AE) of the Gamow- Teller 
strength in this channel. That actually bodes well for any 
attempt to use the nuclear systematics of lighter nuclei 
like 28 Si and 56 Fe to effect an extrapolation of neutrino 
pair emission mechanisms to the higher mass nuclei of 
most interest in stellar collapse. This is because the spin- 
orbit splitting is relatively constant across nuclear mass 
in the range over which we are interested (47j . Particle- 
hole repulsion probably plays a lesser role than spin-orbit 
splitting in pushing strength to larger AE. Interestingly, 
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FIG. 7: 28 Si strength and neutrino pair spectrum at 30 MeV 
excitation computed with spin-orbit splitting in the nuclear 
Hamiltonian set to zero. The strong central peak, lack of 
wings, and concomitant drastic reduction in neutrino pro- 
duction confirm the role of spin-orbit splitting and spin-flip 
transitions in producing the high rates computed in our more 
realistic models. 



Nucleus 


Ji 


Excitation 


Rate 


T (A=56) 


28 Si 


1-4 


10.8 


0.30 


1.2 




0, 2-5 


14.0 


0.02 


1.4 




0-5 


15.0 


0.22 


1.5 




5 


20 


0.50 


1.7 




0-5 


20.0 


0.50 


1.7 




0, 2-5 


21.9 


1.1 


1.8 




5 


25 


1.8 


1.9 




1-5 


25.0 


1.8 


1.9 




0-5 


27.6 


3.6 


2.0 




5 


30 


3.4 


2.1 




0-3, 5 


30.0 


3.3 


2.1 




1-5 


31.5 


4.0 


2.1 




5 


35 


4.5 


2.2 




5 


40 


6.0 


2.4 




0-5 


40.0 


8.6 


2.4 


s Si (no SO) 


5 


30 


0.21 


2.1 


29 Si 


11/2 


30 


4.6 


2.1 


28p 


5 


30 


9.7 


2.1 


4 'Ti 


3/2 


23.0 


0.33* 


1.8 




3/2, 5/2 


27.6 


2.0* 


2.0 


bb Fe 


2 


10.0 


0.01* 


1.2 







15.0 


0.09* 


1.5 




1 


20.0 


2.7* 


1.7 




2 


25.0 


0.30* 


1.9 




0, 1 


27.6 


4.6 


2.0 







30.0 


10.9* 


2.1 




4 


35.0 


12.1* 


2.2 



TABLE I: Energy loss rate (MeV/s/baryon) for various nuclei 
as functions of excitation energy (MeV) and the correspond- 
ing (rescaled to A=56) temperature (MeV). The unstarred 
rates were computed by averaging over several (5 to 14) ini- 
tial states at the indicated excitation energy; starred rates 
were computed from 1 or 2 states. The angular momenta 
of all initial states are indicated. The entry 28 Si (no SO) 
was computed by neglecting spin-orbit splitting in the nuclear 
Hamiltonian. 



the particle-hole repulsion we find in our shell-model cal- 
culations may have a direct analog in the bulk matter 
renormalization of the energetics of weak interaction pro- 
cesses as found in Ref. [481 ] . 

There are three obvious effects of skewing the strength 
distributions to higher transition energy AE. These 
follow from the simple fact that the neutrino-pair en- 
ergy emission rate derived from the strength function is 
weighted by six powers of AE. First, more strength at 
higher AE generally means faster neutrino-pair emission 
rates. In turn, this means more energy will be pumped 
into neutrino pairs by this process. Second, the neutri- 
nos and antineutrinos produced by this process will have 
higher energies on average. This brings up an obvious 
question: will the neutrino pair energies now be so high 
(> 10 MeV) that they are more readily trapped? Third, 
more configurations in play and more configuration mix- 
ing at higher excitation energy will make the transition 
strength and the neutrino-pair emission rate more sensi- 
tive to temperature. 



Higher temperatures and higher excitation energies 
bring up an issue which is unresolved in our work. When 
nucleons are promoted into the next-higher oscillator 
level, how is the Gamow- Teller strength affected? For ex- 
ample, for 28 Si the actual level density just above about 
10 MeV will start to be dominated by one-nucleon excita- 
tions into and out of the sd shell (negative-parity states) . 
The Gamow- Teller down-strength for these states will be 
reduced since they must go to negative-parity final states 
that only exist at a higher energy than exist for posi- 
tive parity. Starting at 20 MeV, more positive parity 
states can be made from two nucleons excited into and 
out of the sd shell. Does this result in a lower overall 
amount of strength? Or is this loss of strength compen- 
sated by more transitions between spin-orbit partners in 
the higher energy shell? In general, it may seem rea- 
sonable that the higher the temperature and excitation 
energy, the "looser" the nucleus and the more transi- 
tions are unblocked [5]. However, this trend could be 
thwarted by the actual behavior of the level density in 
the shell model. Certainly, model space truncation could 
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contribute to this if, for example, not all spin-flip tran- 
sition channels in the higher oscillator level are included 
in the calculation. This issue may or may not compli- 
cate extrapolation of our trends in weak strength energy 
distributions to higher mass nuclei and requires further 
investigation. 

Fig.[5]shows a comparison between our shell- model cal- 
culations of the neutrino-pair energy emission rate per 
baryon for 28 Si as a function of (rescaled) temperature, 
other estimates for this rate in 56 Fe as a function of 
temperature, the neutrino pair emission rate for electron 
bremsstrahlung as a function of temperature assuming 
a density pi2 = 1 [Ijjj], and the neutrino pair emission 
rate for nuclcon bremsstrahlung in nuclear matter as a 
function of temperature [49j]. We include the nucleon 
bremsstrahlung result to show that collectivity within the 
nucleus enhances the emission rate relative to bulk nu- 
clear matter. We found that in the temperature regime 
of interest, our calculations can yield neutrino-pair en- 
ergy emission rates that equal or exceed earlier estimates 
at all temperatures. As a general rule, our estimates 
of these rates are approximately three orders of magni- 
tude faster than neutrino pair production from electron 
bremsstrahlung. 




0.5 1 1.5 2 2.5 



k B T (MeV) 

FIG. 8: Energy per baryon emitted in neutrino pairs as a func- 
tion of temperature. Solid: 28 Si as computed in this paper. 
Dashed: 56 Fe from Fuller & Meyer analytic approximation 
Dotted: 56 Fe from Fuller & Meyer independent single-particle 
shell- model calculation Fine dotted: 56 Fe from Kolb & 
Mazurek Q|. Long dash-dotted: Electron bremsstrahlung into 
neutrino pairs for pi2 = 1 from Dicus et al [23| . Short dash- 
dotted: Nucleon bremsstrahlung into neutri no p airs at nu- 
clear matter density from Friman & Maxwell [491 ] . 

Clearly, de-excitation of nuclei is the dominant contrib- 
utor of relatively low-energy neutrino pairs under these 
conditions. Moreover, the rates presented here are lower 
bounds on the actual neutrino-pair production rates, par- 
ticularly at temperatures between 1 and 1.5 MeV. The 
small number of nucleons in a 28 Si nucleus gives a rela- 
tively low density of states at the temperatures of inter- 
est. As a consequence, there are few lower-lying states to 



transition to, which reduces the total transition rate. In- 
deed, the apparent decrease in emission rate of 28 Si at a 
(rescaled) temperature of 1.4 MeV is a consequence of the 
simple fact that no states near 14 MeV excitation have 
transitions with energies greater than those available to 
states near 10.8 MeV transition. 

// it turns out that the neutrinos produced by pair de- 
excitation have low enough energies to escape the pre- 
supernova star, then this process likely acts as a ther- 
mostat for the collapsing core. In this limit, as the core 
heated up, more neutrino pairs would be produced and 
escape, carrying away entropy, and perhaps keeping the 
core temperature near T — 1 MeV to 1.5 MeV. 

However, our calculations may be suggesting that 
the neutrino pairs produced by de-excitation at higher 
temperature are so energetic that they do not escape. 
Though the thermostat effect will be disabled in this case, 
lower to intermediate energy neutrino phase space will be 
filled more quickly by this process, and the core will ap- 
proach beta equilibrium sooner. This effect would tend 
to block electron capture and neutronization sooner also, 
but against this neutrino-nucleus and neutrino-electron 
down scattering will tend to heat the system, adding en- 
tropy, implying faster electron capture through more free 
protons and nuclear thermal unblocking |5j and, there- 
fore, a lower Y e , a smaller homologous core, and a con- 
comitantly lower initial bounce shock energy. 

The issue of higher energy neutrino-pairs is compli- 
cated further when considering the effects of dynamics 
and flavor. The energy is shared between a neutrino and 
an anti-neutrino, and the energy is small compared to the 
mass of the nucleus. As a consequence, the energy can be 
shared unequally between the two neutrinos, with the nu- 
cleus absorbing whatever momentum is needed to satisfy 
conservation. So a low-energy partner to a high-energy 
neutrino could escape. Furthermore, this process is flavor 
blind: it produces neutrinos of all flavors at equal rates. 
Because high energy neutrinos within the collapsing core 
are produced primarily from electron capture, there is 
plenty of available phase space for nuclear de-excitation 
to produce a high energy electron anti-neutrino or a high 
energy neutrino or anti-neutrino of mu or tau flavor with 
a low energy partner that easily escapes. The resultant 
asymmetry between electron flavored neutrinos and anti- 
neutrinos could impact lepton number within the core, as 
it is more likely to produce a high energy electron anti- 
neutrino with a low energy partner than the other way 
around. 

Another effect of the larger AE values suggested 
by our shell-model calculations may be an enhance- 
ment in the plasmon-mediated neutrino-pair nuclear de- 
excitation process pointed out by Horowitz in Ref. Q. 
The matrix element for the first forbidden vector chan- 
nel considered in Ref. Q is oc (f\q ■ T z \i), where q ~ AE 
is the momentum transfer. This first forbidden channel is 
in general cut down by a geometric factor, (qR) ~ 1/16, 
which is the square of the ratio of the nuclear radius R 
to the inverse momentum transfer. Our larger values of 
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AE should give a smaller reduction, increasing the over- 
all rate of nuclear de-excitation into pairs. 

Though our shell-model calculations are only a begin- 
ning, they do suggest that nuclear de-excitation into neu- 
trino pairs is likely the dominant source of low to interme- 
diate energy neutrino pairs in stellar collapse. Our calcu- 
lations suggest a spin-orbit splitting-induced increase in 
the rate of this process and a steepening of the temper- 
ature dependence of this rate. These calculations also 
suggest, however, that the neutrinos produced in this 
process are more energetic and may be trapped. Only 
inclusion in a full core collapse neutrino transport sim- 
ulation could reveal what role this process plays in core 
collapse supernova explosions. 
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